Evolution of Ultracold, Neutral Plasmas 



O 

o 

> 
O 



c*3 

1/3 

o 

C/5 



> 

On 
^1- 



S. Mazevet, L. A. Collins, J. D. Kress 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

(submitted to PRL) 



CZ3 

o 

>>: 
^ : 

Or 



X 



We present the first large-scale simulations of an ultracold, 
neutral plasma, produced by photoionization of laser-cooled 
xenon atoms, from creation to initial expansion, using classi- 
cal molecular dynamics methods with open boundary condi- 
tions. We reproduce many of the experimental findings such 
as the trapping efficiency of electrons with increased ion num- 
ber, a minimum electron temperature achieved on approach 
to the photoionization threshold, and recombination into Ry- 
dberg states of anomalously-low principal quantum number. 
In addition, many of these effects establish themselves very 
early in the plasma evolution (~ ns) before present experi- 
mental observations begin. 



That a common characteristic connects such diverse 
environments as the surface of a neutron star, the initial 
compression stage of an inertial confinement fusion cap- 
sule, the interaction region of a high-intensity laser with 
atomic clusters, and a very cold, dilute, and partially- 
ionized gas within an atomic trap, seems at first rather 
remarkable. Yet all these cases embrace a regime known 
as a strongly-coupled plasma. For such a plasma, the in- 
teractions among the various constituents dominate the 
thermal motions. The plasma coupling constant r Q , the 
ratio of the average electrostatic potential energy Z a /& a [ 
Z Q , the charge; a Q , the ion-sphere radius = 3/(47rn Q 1 / 3 ); 
and n a , the number density for a given component a] 
to the kinetic energy ksT Q , provides an intrinsic mea- 
sure of this effect [[!]. When r Q exceeds unity, vari- 
ous strong-coupling effects commence such as collective 
modes and phase transitions. For multi-component plas- 
mas, the coupling constants need not be equal or even 
comparable, leading to a medium that may contain both 
strongly- and weakly- coupled constituents. 

Since temperatures usually start around a few hundred 
Kelvin, most plasmas found in nature or engineered in 
the laboratory attain strongly-coupled status from high 
densities, as in the case of a planetary interior or a shock- 
compressed fluid or fr° m highly-charged states, as 
in colloidal or "dusty" plasmas B. In both situations, 
the particle density usually rises well above 10 18 /cm 3 . 
On the other hand, ion trapping and cooling methods 
have produced dilute, strongly-coupled plasmas by rad- 
ically lowering the temperature. At first, these efforts 
were limited to nonncutral plasmas confined by a mag- 
netic field ||. However, recently, new techniques 
have generated neutral, ultracold plasmas, free of any 



external fields at densities of the order of 10 8 -10 10 /cm 3 
and temperatures at microkelvins (^K). 

Two methods, one direct and one indirect, but both 
employing laser excitation of a highly-cooled sample of 
neutral atoms, have successfully created such neutral 
plasmas. The direct approach employs the photoioniza- 
tion of laser-cooled xenon atoms ^)-^| , while the indirect 
generates a cold Rydberg gas in which a few collisions 
with room-temperature atoms produces the ionization 
P,M- In both cases, the electron and ion temperatures 
start very far apart with the former from IK to 1000K 
and the latter remaining at the initial atom temperature 
of a few /iK. This implies that the coupling constant for 
the electrons (r e ) ranges between 1 and 10, while the 
ions begin with Tj ~ 1000. Based on the results pre- 
sented below, the following picture emerges. As the sys- 
tem evolves, the electrons reach a quasi-equilibrium while 
at the same time beginning to move the ions through 
long-range Coulomb interactions. Eventually, the ions 
also heat, and the whole cloud begins to expand with 
the ions and electrons approaching comparable temper- 
atures. The former processes occur on the order of pico- 
to nano- seconds, while the full expansion becomes no- 
ticeable on a microsecond scale. Therefore, by following 
the progress of the system, we can study not only the 
basic properties of a strongly-coupled plasma but also its 
evolution through a variety of stages. The path to equi- 
libration still remains a poorly understood process for 
strongly-coupled plasmas in general. 

The dilute character of the ultracold plasmas pro- 
vides a unique opportunity to explore the intricate nature 
of atomic-scale, strongly-coupled systems, both through 
real-time imaging of the sample to study wave phenom- 
ena, collective modes, and even phase transitions, and 
through laser probes to examine internal processes. The 
initial experiments have already discovered unexpected 
phenomena such as more rapid than expected ion ex- 
pansion 0] and Rydberg populations that strongly de- 
viate from those predicted from standard electron-ion 
recombination rates |q JQ]-|l3|] . So far, the experiments 
have examined the later expansion stage of the whole 
plasma cloud, which occurs on the scale of microseconds. 
However, the nature of the plasma at this stage depends 
strongly on its evolution from its creation. While ultra- 
fast laser techniques offer the eventual prospect of prob- 
ing these early times, for now, only simulations can pro- 
vide an understanding of the full evolution of these sys- 
tems. We focus on the first such generated plasma |J 
from the photoionization of cold Xenon atoms; however, 



the findings also provide insight into other ultracold sys- 
tems. 

Since the electron temperature greatly exceeds the 
Fermi temperature for this system, we may effectively 
employ a two-component classical molecular dynamics 
(MD) approach in which the electrons and ions are 
treated as distinguishable particles, interacting through a 
modified Coulomb potential. We have employed several 
different forms from an ab-initio Xe + + e~ pseudopoten- 
tial [jlij to a simple effective cut-off potential based 
on the de Broglie thermal wavelengths; we found little 
sensitivity of the basic properties to this choice. This 
finding stems from the dominance of the interactions by 
the very long-range Coulomb tails since the average par- 
ticle separations are of the order of 10 4 A(1 to 2 times the 
de Broglie wavelength for electrons at a temperature of 
O.IK). The short-range modification of the Coulomb po- 
tential basically prevents particles with opposite charges 
from collapsing into a singularity. Due to the open and 
nonperiodic boundary conditions, sample size effects can 
play a critical role in the simulations Jl6| ]. Therefore, we 
must employ particle numbers of the basic order of the 
experiments, requiring efficient procedures for handling 
the long-range forces. In addition, the extended simula- 
tion times needed to model the entire plasma evolution 
[fs to /is] demand effective temporal treatments. To this 
end, we have employed multipole-tree procedures [|l7j in 
conjunction with reversible reference system propagator 
algorithms p) (r-RESPA) through the parallel 3D MD 



program NAMD |19|, generally used in biophysical ap- 
plications. These procedures allow us to treat samples of 
between 10 2 and 10 4 particles [N = N 4 /2 = N e /2]; most 
simulations were performed with N = 10 3 . 

To reproduce the initial experimental conditions || , we 
distribute the ions randomly according to a 3D Gaussian, 
whose rms radius, tr, matches the desired ion number 
density ni [ni = Ni/(27rcr 2 ) 3 / 2 ] and impose a Maxwell- 
Boltzmann velocity distribution at Ti — 1/iK. To each 
ion, we associate an electron placed randomly in an orbit 
of radius r = 50A. The electron velocities point in ran- 
dom directions but with fixed magnitude, determined by 
the photoionization condition [K* = ^m e w 2 — Ticu — l/r , 
with u>, the laser frequency and v\ as a sum over all three 
spatial components]. This idealized prescription models 
the photoionization process by allowing each electron to 
escape the Coulomb well with a final prescribed energy. 
By varying the laser frequency, we can control the final 
effective electron temperature, just as in the experiments. 
We tested this particular ionization model against simu- 
lations where both the initial electron-ion radius r„ and 
the form of the ion distribution are varied. Overall, we 
find little sensitivity to these variations on the basic ini- 
tial conditions. Finally, we take a kinematic definition 
of temperature as the average kinetic energy per parti- 



The plasma evolves through several stages. In the first 
or photoionization stage, which lasts on the order of fs, 
the electrons climb out of the Coulomb well and become 
basically freely-moving particles. In the next stage, the 
electrons reach a quasi-equilibrium at T e due to their 
fast intra-particle collisions. Following this stage, the 
electron-ion collisions begin the slow process of heating 
the ions, which can require up to /xs. This process can 
also be viewed in terms of an electron pressure term in a 
hydrodynamical formulation |t]]. Then, the whole cloud 
of ions and electrons begins a systematic expansion. This 
progression clearly evinces itself in Fig. 1, which shows 
the temporal behavior of electron and ion average kinetic 
energies, using an effective mass for the ion of m; = 0.01 
amu in order to accelerate the evolutionary process. The 
averaged kinetic energy of the electrons and ion become 
comparable and the cloud perceptibly expands in about 



20 ns. Scaling this value by y' yields an estimate for 

this time in a Xe plasma of about l/is, in line with the 
experiments ||. Also consistent with the experimental 
measurements |Q , these preliminary calculations indicate 
that a large fraction of the kinetic energy transferred by 
the electrons results in outward translational motion for 
the ions. This implies that the kinetic temperature de- 
fined above coincides with the usual themodynamic tem- 
perature (random thermal motion) for the electrons only. 
In all other simulations discussed, the ions carry the mass 
of Xe. 
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FIG. 1. Schematic variation of the electron (open square) 
and ion (filled square) average kinetic energy per particle as 
a function of time for a model plasma with m;=0.01 amu. 



Since we shall use the simulations to understand in- 
termediate stages in the development of the ultracold 
plasma, we need to establish the validity of the procedure 
by comparing to experiments. Two initial observations 
have particular significance: 1) the number of trapped 
electrons rises with the number of ions created for a fixed 
final temperature, and 2) the electrons attain a certain 
minimum temperature no matter how small Tilo becomes. 
Figure 2a shows that the number of trapped electrons 
as a function of the number of ions and initial electron 



temperature basically follows the general experimental 
trends Bj. Some of the electrons, freed by the photoion- 
ization process from ions with sufficient energy, escape 
the atomic cloud and never return, leaving an overall 
positively charged system. This residual charge then ef- 
fectively traps the remaining electrons so that the cen- 
ter of the distribution resembles a neutral plasma. The 
larger the number of ions produced, the more effective 
the confinement of the electrons. 

In Fig. 2b, we present for nj=4.32xl0 9 ions/cm 3 , the 
electron temperature T{ after the initial equilibration in 
the simulations as a function of the excess initial temper- 
ature T*[= g-K"e] given to the electrons during the pho- 
toionization process. The plateau at about 5K [T e ~ 1] 
for small T* appears quite pronounced. This effect arises 
due to a mechanism usually designated as "continuum 
lowering" [EI. Despite the small energies and enormous 
distances involved, at least by usual atomic standards, 
the interaction of an ion with its neighbors still has a no- 
ticeable effect due to the very long-range of the Coulomb 
potential. This shifts the appropriate zero of energy from 
the isolated atom to the whole system. Using a simple 
binary interaction, we can estimate this energy differ- 
ence for an electron halfway between the ions [~ 2/ai on ] 
at around 3K for this density, in qualitative agreement 
with the simulations. Even for a photon energy near the 
atomic threshold, the electron within the plasma will still 
gain this minimum energy. This same effect also explains 
the dependence of the final electron temperature on the 
ion density. In this situation, an increase in the ion den- 
sity enhances the zero of energy shift and leads to an 
increase in the final electron temperature. 
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FIG. 2. (a) Percentage of electrons trapped as a function of 
the number of ions created for two electron temperatures (3K 
and 10K). (b) Variation of the final electron temperature T{ 
as a function of the initial excess kinetic energy represented 
an associated temperature T*. 




experiments H have used this quantity as an indirect 
means of deducing the density of the ultracold plasma. 
Given the nature of the system with open boundary con- 
ditions and quasi-equilibration among the electrons, the 
question arises as to whether f e has a precise definition. 
Neglecting temperature dependence, the electron plasma 
frequency is given by 



(1) 



where e is the elementary charge, n e the electron den- 
sity, and m e the electron mass. From the MD simula- 
tions, the electron plasma frequency is obtained from the 
Fourier transform of the electron velocity autocorrela- 
tion function pofl . We find, reassuringly, that a distinct 
though broadened peak arises near the f e predicted by 
Eq.|l| in the regime in which the electrons reach a quasi- 
equilibration. For example, at a density of 4.32xf0 10 
ions/cm 3 , an electron temperature of 3K, and simula- 
tion time extending up to 18ns, MD gives plasma fre- 
quencies of 1.2GHz (1.4 GHz) while Eqjl] yields 1.57GHz 
(1.87GHz) for trapped electrons at 70% (100%) of the ion 
density. When the ion density is decreased to a value of 
4.32xl0 9 ions/cm 3 at the same electron temperature, / e 
from the MD simulations becomes 0.5Ghz, also in accor- 
dance with Eq.[l] . In general, we find that the number of 
particles used in the simulation cell has little effect on the 
determination of the electron plasma frequency beyond 
providing a better statistical sample. 
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FIG. 3. Time evolution of the modified ion-electron pair 
correlation function Gi e as a function of density and electron 
temperature: (a)-(c) 4.32xl0 10 ions/cm 3 and 10K, and (d)-(f) 
4.32xl0 9 ions/cm 3 and 3K. 



The molecular dynamics simulations yield additional 
properties of the ultracold plasma. We paid particu- 
lar attention to the electron plasma frequency f e since 



We now turn to the nature of the constituents as the 
plasma evolves and pay particular attention to the for- 
mation of Rydberg atoms as noted in the experiments 



|^| jTo|. To this end, we examine over space and time 
the ion-electron pair correlation function gi e (r), which 
gives the probability of finding an electron a distance r 
from a reference ion p"5| . To enhance the examination of 
the closer encounters, we have employed a modified form 
Gi e (r) by subtracting a uniform distribution and aver- 
aging <7i e (r) over a small time interval of 1.5 ns. Figure 
3 displays Gi e (r) for two representative cases. The open 
boundary conditions, which force the distribution to zero 
at some large distance, give Gi e (i) a distinctly different 
behavior than seen in a fluid or periodic system. 

Figures 3(a)-(c) display the development of a plasma 
for ni=4.32xl0 10 ions/cm 3 and T e ~ 10K in the regime in 
which the electrons have reached quasi-equilibration and 
the ions have started to heat. After 2.5ns, the electron 
distribution departs from its initial form, especially be- 
yond a radius of 10000A. By 8.5ns, we begin to observe an 
increase in Gi e just below 1000A. This signifies the forma- 
tion of Rydberg states in an average principal quantum 
number n ~ 40, a conclusion supported by an examina- 
tion of movies of the simulation. These movies also show 
stable Rydberg atom configurations over many orbital 
periods of the electrons. A simple model (l3) that bal- 
ances various collisional processes, using strictly atomic 
cross sections, predicts a much larger value (n > 100). 
For conditions closer to the experimental measurements 
(lower density), as depicted in Figs. 3(d)-(f), we find a 
significant Rydberg population after only 20ns (Fig.3f). 
The electron distribution shifts to around 2000A, corre- 
sponding to a principal quantum number n ~ 60. In both 
cases, the Rydberg population is estimated to be between 
5% and 10 % of the total number of electrons. These 
findings closely resemble those of the experiments [||, 
though for later times (~ /is) and indicate that collective 
or density effects may play an important role by changing 
the accessible Rydberg-level distribution or cross sections 
through long-range interactions. By examining the elec- 
tron mean square displacement over time, we can identify 
two stages in the evolution of the cold ionized gas. First, 
after the ionization of the atoms, the electrons diffuse 
throughout the system for several nanoseconds with no 
significant Rydberg atom formation. Second, the elec- 
trons reach the edge of the cloud and begin systematic 
multiple traverses of the system. The Rydberg atom pop- 
ulation becomes noticeable only several nanoseconds into 
this second stage. The feature that attracts the most in- 
terest is the relatively short time scale required to estab- 
lish the neutral plasma and produce a noticeable Rydberg 
population. 

In summary, we have followed the evolution of an ul- 
tracold, neutral plasma over a broad range of tempo- 
ral stages with classical molecular dynamics simulation 
techniques. We find general agreement with experimen- 
tal observations of the number of trapped electrons, the 
minimum of electron temperature, and the production of 
Rydberg atoms in low-lying states. The latter two con- 



ditions especially demonstrate the importance of strong- 
coupling or density effects on the basic atomic inter- 
actions. In addition, we have found that the electron 
plasma frequency appears as a valid tool to probe the 
state of the system. Important to the understanding of 
the temporal development of these plasmas, we discov- 
ered that recombination and the formation of long-lived 
Rydberg states occur rapidly on the order of nanosec- 
onds. Our studies continue in an effort to push into the 
fully expanded regime by using hydrodynamical methods 
tied to the molecular dynamics. 
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